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iSTRACT 


Analyses  of  trace  and  major  elements  in  gronp  III-V  compound 


matrices  by  secondary  ion  mass  spectrometry  (SIMS)  have  shown  that 


practical  ion  yields  are  1 'nearly  related  to  the  matrix  composition. 


The  affinity  of  the  matrix  to  oxygen  appears  to  be  the  critical  factor 


in  this  relationship.  A  direct  relationship  between  the  slopes  of 


these  calibration  lines*  determined  for  elements  in  the  same  column  of 


the  periodic  table*  and  the  first  ionization  potential  of  the 


respective  elements  has  also  been  shown. 
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Secondary  ion  aaes  spectrometry  (SIMS)  is  a  powerful  technique 
for  the  analysis  of  thin  solid  films.  It  has  ppm  sensitivity  for  most 
elements,  excellent  depth  resolution  (<100  A),  and  good  lateral 
resolution  (about  1.0  pm).  Using  standards  prepared  by  ion 
implantation,  quantitative  analysea  accurate  to  about  15%  can  be 
obtained  for  trace  elements  (<0.1  atomic  percent)  in  homogeneous 
matrices  [1-3].  Due  to  matrix  effecta.  however,  the  standard  must 
match  the  sample  matrix.  This  not  always  convenient  or  possible.  In 
addition.  the  quantification  of  SIMS  depth  profiles  in 
multilayer-multimatrix  samples  remains  a  problem  due  to  the  changing 
practical  ion  yields  and  sputtering  yields  which  are  encountered. 

For  Al^a^xAa  and  related  matrices,  it  has  recently  been  shown 
that  relative  ion  yields  (Ex)  and  relative  sputtering  yields  (MS)  vary 
linearly  with  sample  composition  [4] .  The  relative  values  were 
obtained  by  normalising  practical  ion  yields  (r)  and  sputtering  yields 
(S)  from  a  sample  matrix  to  those  from  a  standard  matrix  when  both 
measurements  were  performed  under  nearly  identical  analysis 
conditions.  By  this  procedure,  precise  calibration  lines  have  been 
obtained,  and  subaequently  been  applied  to  a  variety  of  AlzGaj_zAs 
superlattices  using  the  depth  profile  correction  program  SLIC 
(superlattiee  and  interface  calibration)  [5]. 

The  general  application  of  such  a  calibration  method  to  a  greater 
variety  of  matrices  requires  an  understanding  of  the  fundamental 
processes  involved.  There  are  currently  two  principal  hypotheses  for 
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the  variation  of  t'l  with  matrix.  It  has  been  suggested  that  matrix 
effects  are  primarily  a  function  of  sputtering  yields  [6,  7].  Lover 
S' a  enhance  the  build-up  of  reactive  primary  ions  (0|  or  Cs+)  in  the 
surface  resulting  in  increased  t's.  This  hypothesis  is  expressed  in 
eq.  1 


Etg  «c  <l/RS)y 


(1) 


vhere  y  is  a  constant  determined  experimentally  for  the  analyte  g. 
Others  have  asserted  that  under  certain  conditions  r's  are  a  linear 
function  of  matrix  composition  [8,  9].  This  compositional  approach  is 
expressed  in  eq.  2 


■'i  -  <2> 

vhere  n  is  the  number  of  elements  in  the  matrix,  is  the  atomic 
fraction  of  element  i  in  the  matrix,  and  P^g  is  a  dimensionless 
parameter  representing  the  influence  of  element  i  on  the  ion  yield  of 
element  g.  In  addition,  a  strong  correlation  betveen  the  mean  free 
energies  of  matrix-oxygen  bonds  and  the  observed  trends  in  ionization 
probabilities  has  been  reported  for  certain  binary  metal  alloys  [10, 
111. 
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la  the  present  study,  these  principal  explanations  for  the 
variation  of  r's  with  natrix  were  investigated  nsing  Rr's  from  trace 
and  major  elements  in  various  group  I1I-V  compound  matrices.  A  strong 
relationship  between  R-c's  and  RS's  was  not  observed.  Alternatively,  a 
very  strong  linear  relationship  was  found  between  the  Rr's  and  the 
average  bond  energies  of  the  saaple  matrices  to  oxygen.  Similar 
linearity  was  found  during  Ar+  bombardment  with  a  high  ambient  oxygen 
pressure  in  the  sample  chamber.  For  elements  in  the  same  column  of 
the  periodic  table,  a  direot  correlation  was  observed  between  the 
slopes  of  these  lines  and  the  ionization  potentials  of  the 
corresponding  analytes.  Using  these  trends,  it  is  possible  to  predict 
when  matrix  effects  will  be  a  problem.  Then,  using  the  appropriate 
calibration  lines  it  is  possible  to  correct  for  these  matrix  effects. 


Sasrole  Preparation.  The  Al^Ga^jAs,  Ga  ^In  ^As,  GaSb,  and  InSb 
matrices  were  grown  by  molecular  beam  epitaxy  (MBE)  on  semi-insulating 
GaAs,  InP,  GaAs,  and  InSb  substrates,  respectively.  Semi-insulating 
InP  and  GaP  substrates  were  also  used  as  sample  matrices.  The  matrix 
compositions  were  determined  from  MBE  growth  parameters  and  verified 
to  an  accuracy  of  better  than  10%  [12]  using  Rutherford  backscsttering 
(RBS)  and  photoluminescence  spectroscopies.  Trace  elements  were 
introduced  using  ion  implantation.  Generally,  two  different  elements 
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wore  iaplanted  into  each  sample.  Prior  to  implantation  samples  were 
cleaned  with  trichloroethylene.  The  implantation  parameters  are 
listed  in  table  1. 

Instrnmentation.  SINS  analysis  was  carried  out  on  a  CAMECA  (103. 
Boulevard  Saint-Denis.  92403  Courbevoie  Cedez.  France)  IMS-3F  ion 
microanalyzer  using  an  electron  multiplier  in  the  pulse  counting  mode 
for  signal  detection  [13].  The  instrument  comes  interfaced  to  a 
HEWLETT  PACKARD  (3404  East  Harmony  Road.  Fort  Collins.  Colorado  80525) 
9845T  microcomputer  for  control  and  data  acquisition.  The 
experimental  parameters  are  listed  in  Table  2.  To  avoid  saturating 
the  electron  multiplier*  the  sampling  area  was  reduced  to  5.0  X  10~7 
cm^  when  ^Al+,  *®Ga+,  ^Li+,  and  ^*Mg+  were  analyzed.  A  multiple 
sample  holder  was  used  to  simultaneously  mount  several  samples.  Depth 
measurements  on  the  sputtered  craters  were  performed  using  a 
mechanical  stylus. 

Software.  Programs  for  instrumental  control,  data  analysis,  and 
matrix  correction  were  written  in  BASIC  for  the  HEWLETT  PACKARD  9845T. 

Procedure.  Following  ion  implantation,  4  to  7  samples  were 
simultaneously  mounted  and  depth  profiled  using  a  multiple  sample 
holder.  Each  group  of  samples  included  GaAs  as  a  standard  matrix  for 
ion  yield  and  sputtering  yield  normslization.  The  ssmples  were 
inserted  simultaneously  to  insure  nearly  identical  analysis 
conditions.  After  allowing  the  pressure  in  the  sample  chamber  to 
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resell  a  steady  state  condition,  the  02+  primary  ioa  beam  was  focused, 
sad  the  proper  mass  settiags  were  determined.  Without  maaipulatiag 
say  iastruaeatal  parameters,  the  samples  were  aaalyzed  eoaseeutively 
uatil  at  least  3  profiles  of  each  sample  had  beea  made,  x's,  Ex's, 
S's,  sad  ES's  were  determiaed  usiag  staadard  procedures  [4].  The  x's 
sad  Ex's  for  ^B,  ^lp#  sad  ^As  were  also  determiaed  usiag  aa  Ar+ 
primary  ioa  beam  while  the  sample  chamber  was  flooded  with  oxygea  to  a 
pressure  of  3  X  10~*  torr. 

RESULTS  AND  DISCUSSION 

Siace  Ex's  are  liaearly  related  to  both  ES's  sad  sample 
eompositioa  for  Al^a^^As  matrices  14],  the  alteratioa  of  x's  with 
matrix  caa  be  explaiaed  by  either  of  the  two  hypotheses  of  siatrix 
effects.  The  validity  aad  applicability  of  these  two  caa  only  be 
evaluated  by  exteading  the  calibratioa  to  other  group  III-V  compound 

matrices.  The  sputtering  yield  hypothesis  was  evaluated  directly. 

/ 

The  compoaitioaal  approach  was  evaluated  usiag  the  observation  of  Tu 
aad  Eeuter  [10,  11]  that  the  affinity  of  a  matrix  for  oxygen  (AG£ 

metal  oxide)  determines  the  influeace  of  a  particular  matrix 
composition  on  x's.  Consequently,  the  parameter  P^  g  was  treated  as  a 
measure  of  the  affinity  of  element  i  for  oxygen. 

The  exact  forms  of  the  metal-oxygen  complexes  which  are  created 
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during  the  sputtering  process  ere  difficult  to  determine. 
Consequently*  the  free  energies  of  metal  oxide  formation  were  not  used 
as  an  indicator  of  oxygen  affinity.  Instead,  this  affinity  was 
roughly  approximated  using  the  bond  energies  of  the  diatomic  metal 
oxides.  Preliminary  analysis  indicated  that  a  linear  relationship 
existed  between  x'a  and  the  oxygen  affinity  of  the  matrices,  but  that 
this  linearity  could  be  improved  for  all  analytes  by  using  slightly 
modified  bond  energy  values  as  an  indicator  of  oxygen  affinity.  These 
modified  bond  energy  values  are  used  throughout  the  paper.  The 
literature  values  [14]  and  the  modified  values  of  these  bond  energies 
are  given  in  Table  3. 

The  two  approaches  to  the  problem  of  matrix  effects  are  compared 
in  Figures  1-4  for  28Si,  9Be,  31P,  and  24Mg  in  various  group  III-V 
semiconductors.  According  to  the  compositional  theory,  plots  of  Sr 
versus  the  average  matrix-oxygen  bond  energy  should  yield  straight 
lines.  Similarly,  according  to  the  sputtering  yield  hypothesis,  plots 
of  log(Rx)  versus  log(l/KS)  should  also  yield  straight  lines.  As 
apparent  in  these  figures,  the  compositional  theory  satisfies  this 
criterion  much  better  than  the  sputtering  yield  approach.  In  fact, 
the  linearity  that  is  apparent  in  the  plots  of  log (Ex)  versus 
log(l/RS)  can  be  attributed  almost  entirely  to  the  linearity  expected 
for  the  Al^aj-xAs  data  points.  The  linear  correlations  and  relative 
standard  deviations  of  the  slopes  of  these  two  types  of  plots  are 
preaented  in  Table  4.  The  linear  correlations  of  the  compositional 
approach  are  quite  superior  to  those  of  the  sputtering  yield  approach. 


In  addition,  the  relative  standard  deviations  of  the  slopes  obtained 
by  the  sputtering  yield  aethod  vere  aany  times  larger  than  those 
obtained  by  the  coa positional  aethod. 

Certainly,  these  experimental  data  favor  the  compositional 
approaoh.  If,  however,  aatrix  effects  in  these  saaples  are  determined 
by  the  affinity  of  the  saaples  for  oxygen,  one  would  also  expect  a 
linear  relationship  between  Rt's  and  the  average  matrix-oxygen  bond 
energies  when  an  Ar+  primary  ion  beaa  is  used  in  conjunction  with 
oxygen  flooding.  Suoh  an  investigation  was  carried  out  on  ^B,  31p# 
and  ?5As  under  conditions  short  of  oxygen  saturation.  As  apparent  in 
Table  5,  tl  ,  linear  relationship  between  Re's  and  these  bond  energies 
is  aaintained  as  expected  despite  the  different  analysis  conditions. 
In  fact,  as  shown  in  Table  6,  the  slopes  of  the  ^B,  ^lp,  an(j  75^s 
calibration  lines  under  both  types  of  bombardment  are  very  similar. 
The  differences  that  are  apparent  can  most  likely  be  attributed  to  the 
under  saturation  of  the  sample  matrices  with  oxygen.  Vith  the 
addition  of  these  corroborating  data,  it  appears  that  matrix  effects 
in  group  III-V  semiconductors  are  a  linear  function  of  the  sample 
composition  and  in  particular  the  affinity  of  a  particular  matrix  to 
oxygen.  Oxygen  flooding  and  the  0£  bombardment  process  control  the 
availability  of  oxygen,  but  the  oxygen  affinity  of  the  matrix 
determines  the  extent  of  oxidation.  For  Al^a^.jAs  [4,  5,  15]  and 
other  group  III-V  semiconductors,  the  oxygen  affinity  is  the  most 
critical  factor  influencing  matrix  effects. 
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In  addition  to  determining  a  general  mechanism  for  matrix 
effects,  it  is  also  important  to  understand  how  matrix  effects 
inflnenoe  the  t’»  of  one  particnlar  element  versus  another.  One  can 
gain  an  insight  into  this  phenomenon  by  comparing  the  slopes  of  the 
calibration  lines  io  the  first  ionization  potentials  of  the  respective 
analytes.  Such  a  comparison  is  presented  in  Table  6.  At  first 
glance,  there  is  no  apparent  relationship  between  the  slopes  and  the 
first  ionization  potentials.  However,  upon  examining  just  those 
elements  which  fall  in  the  same  column  of  the  periodic  table  [(B,  Al, 
Ga),(P,  As,  Sb),  and  (Be,  Mg)],  a  general  relationship  is  observed. 
Those  elements  with  higher  first  ionization  potentials  yield 
calibration  lines  of  steeper  slope  than  those  with  lower  first 
ionization  potentials.  This  indicates  that  elements  which  have  a 
small  ionization  probability  are  influenced  by  matrix  effects  (the 
amount  of  bound  oxygen)  to  a  greater  extent  than  those  which  have  a 
large  ionization  probability.  As  one  might  expect,  this  relationship 
is  very  similar  to  that  commonly  observed  for  elemental  r's  in  a 
single  matrix  enhance  1  by  oxygen  bombardment  or  flooding:  elements 
that  have  a  small  ionization  probability  are  enhanced  to  a  greater 
extent  by  the  presence  of  oxygen  than  those  that  have  a  large 
ionization  probability.  Thus,  the  relationship  between  the  slopes  of 
these  calibration  lines  and  the  first  ionization  potentials  is 
intuitive. 


Application.  The  excellent  linearity  of  the  Be  versus 
latrix-oxygen  bond  energy  calibration  lines  and  their  relationship  to 
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elemental  ionisation  potentials  can  be  used  to  improve  the  quality  of 
both  qualitative  and  quantitative  SIMS  analysis.  Using  the 
relationships  expressed  in  these  matrix  calibration  lines,  one  can 
anticipate  when  matrix  effects  will  be  a  problem,  and  hov  they  may 
distort  depth  profiles  of  layered  multimatrix  samples.  For  example,  a 
largely  distorted  gaussian  depth  distribution  would  be  expected  for  a 
implant  through  a  GaAs  layer  into  an  Al^aj^As  layer  because 
is  very  sensitive  to  the  presence  of  oxygen,  i.e.  matrix  effects,  and 

^.3G*.7A*  has  a  much  greater  oxygen  affinity  than  GaAs. 
Alternatively,  a  smaller  distortion  would  be  expected  for  an  2*Mg 
implant  and  an  even  smaller  distortion  for  an  7Li  implant  into  such  a 
structure.  The  sensitivity  of  these  elements  to  matrix  effects 
decreases  going  from  11B  to  24Mg  to  7Li.  The  SIMS  analyses  of  these 
implants  are  shown  in  Figure  5. 

When  a  quantitative  SINS  analysis  in  either  a  single  matrix  or  a 
layered  multimatrix  sample  is  desired,  these  calibration  lines  can  be 
quite  valuable.  Since  these  types  of  calibration  lines  are 
reproducible  [4],  dopant  distributions  in  a  homogeneous  group  III-V 
compound  matrix  can  be  quantified  using  the  calibrstion  lines  and  a 
single  standard  prepared  from  any  of  the  group  III-V  compound 
matrices.  There  is  no  need  to  make  a  separate  standard  for  each 
matrix.  These  calibration  lines  can  also  make  the  quantitative  SIMS 
analysis  of  layered  multimatrix  samples  possible.  For  very  complex 
layered  multimatrix  samples,  these  matrix  calibration  lines  may  be 
inserted  into  a  variation  of  the  program  SLIC  [5]  for  matrix 
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correction.  For  simpler  structures  with  just  s  few  well  defined 
interfaces,  one  can  treat  the  interfaces  as  linear  concentration 
gradients  from  one  matrix  to  another.  Both  approaches  are  presented 
in  Figure  6.  A  ®Bo  concentration  plateau  approximately  0.1  pm  wide 
and  1  -  2  X  10*®  atom/cm*  high  was  grown  by  MBE  while  the  matrix  was 
linearly  changed  from  GaAs  to  Al^Ga^As.  In  the  uncorrected  profile, 
the  ®Be  distribution  (dashed  line)  resembles  a  sharp  spike  rather  than 
a  plateau.  In  addition,  the  thickness  and  peak  concentration  of  the 
plateau  can  not  be  determined.  The  shape  of  the  plateau  differs 
slightly  between  the  two  versions  of  the  corrected  profile.  However, 
both  types  of  corrections  allow  a  dramatic  improvement  over  the 
uncorrected  profile  in  the  determination  of  plateau  thickness  and  peak 
concentration. 

In  summary,  the  influence  of  matrix  effects  on  ion  yields  in 
group  III-V  compound  matrices  can  be  precisely  calibrated.  A 
compositional  approach  to  matrix  calibration  based  on  the  oxygen 
affinity  of  the  sample  fits  the  experimental  data  quite  well  while  the 
sputtering  yield  approach  did  not.  The  influence  of  matrix  effects  on 
elements  in  each  column  of  the  periodic  table  was  shown  to  be  directly 
related  to  the  first  ionization  potential  of  the  analytes.  Finally, 
the  relationships  observed  can  be  used  to  improve  the  quality  of  both 
qualitative  and  quantitative  SIMS  analyses. 
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Table  1 

Ion  Iaplantation  Parameters 


Implant 

Element 

Flnenoe 
(atom/ cm2) 

Energy 

(keV) 

Source 

9Be 

1  X  1014 

250 

Be  solid 

1  X  1014 

250 

BP3  gas 

28Si 

1  X  1015 

250 

SiF^  gas 

31p 

1  X  1015 

300 

PF3  gas 

121Sb 

2  X  1014 

250 

Sb  solid 

24Hg 

1  X  1015 

300 

Ng  solid 

7Li 

1  X  1014 

150 

Li  solid 

ky 
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Table  2 

SIMS  Experimental  Parameters 

Primary  Ion:  Oj*,  Ar+ 

Primary  Ion  Energy:  S.5  keV 

Primary  Ion  Current  Density:  -2.5  X  10-2  A/ cm2 
Easter:  300  X  300  jim  or  400  X  400  |tm 
Sampling  Area:  2.83  X  10“ 5  cm2 
Energy  Vindow:  130  eV 

Residual  Sample  Chamber  Pressure:  3  X  10-8  torr 
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Table  3 

Matrix  Eleaent-Oxygen  Bond  Energies 


Matrix 

Eleaent 

Literature 

Valnes 

(kcal/aole) 

Modified 

Valnes 

(kcal/aole) 

A1 

116  A  5 

116.0 

Ga 

68  *  15 

68.0 

In 

in 

67.0 

P 

120  t  4 

115.0 

As 

115  A  3 

115.0 

Sb 

89  A  20 

114.6 

fJIW  TO  TOOT  •.' 


18 
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hr 


Table  4 


Linearity  of  tbe  Spattering  Yield  and 
Coapositional  Approaches  to  Matrix  Calibration 
for  Group  III-V  Con pound  Matrices 
Under  (>2+  Bonbardaent 


Log  (Mr)  Versus  Log  (1/RS) 


Analyte 

Linear 

Correlation 

BSD 

Slope 

% 

2iSi 

0.633 

64.9 

9  Be 

0.713 

54.1 

31p 

0.731 

30.8 

24*s 

0.748 

21.1 

Mr 

Versus  Matrix-Oxygen  Bond  Energy 

Analyte 

Linear 

Correlation 

r* 

RSD 

Slope 

% 

2*Si 

0.999 

0.28 

9Be 

0.998 

0.42 

31p 

0.991 

0.24 

24Mg 

0.993 

0.60 
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Tabic  3 

Linearity  of  Rt  Versa*  Matrix-Oxygen  Bond  Energy 
For  Group  III  and  V  Compound  Matrices 
Under  Ar+  Boabardaent  Vitb  Oxygen  Flooding 


Linear 

Correlation 

ri 


BSD 

Slope 

% 


Analyte 
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Table  6 

Dependence  of  Rt  Versus  Bond  Energy  Line  Slopes 
On  First  Ionization  Potential 


Under  0j+  Boabardaent 


Analyte 

Slope 

11b 

1.90 

11*1 

0.19 

69Ga 

0.25 

3h 

0.93 

75  As 

0.30 

l^Sb 

0.30 

9Be 

1.86 

24Mg 

0.52 

28Si 

1.88 

7Li 

0.10 

Under  Ar+ 

Boabardaent 

1st 

Ionization 


Intercept 

Potential 

(eV) 

-173 

8.30 

-  17 

5.98 

-  22 

6.00 

-  84 

11.02 

-  26 

9.81 

-  26 

8.64 

-169 

9.32 

-  47 

7.64 

-170 

8.15 

-  9 

Oxygen  Flooding 


5.39 
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Figure  5 . 


Figure  6. 


SIMS  analysis  of  ^B,  ^Mg,  and  7Li  inplants 
through  an  GaAs/Al^aj-^s  structure,  (a)  depth 

profiles  of  Hb+  ( - )  and  75As+  ( _ )  {interface 

at  31  tine  units];  (b)  depth  profiles  of  ^*Mg+  ( - ), 

^Li+  ( - ),  and  73As+  (....)  {interface  at  20 

tine  units] . 

SIMS  analysis  of  a  Be  plateau  at  an 

GaAs/Al' jGa( .7)As  interface,  (a)  the  hypothetical 

Be  ( - )  {1  -  2  X  10*8  aton/cn3  peak  cone.} 

and  A1  ("■■ . )  {6.7  X  10^1  aton/cn3  peak  cone.) 

structure;  (b)  uncorrected  profile  of  ®Be+  C - )  and 

75a#  ( - );  (c)  corrected  profiles  of  Be  ( - ) 

{2  X  1018  aton/cn3  full  scale]  and  A1  ( - ) 

{1  X  10«  aton/cn3  full  scale]  obtained  assnning  a 

linear  concentration  gradient  at  the  interface;  (d)  corrected  j 

profiles  of  Be  ( - )  {2  X  1018  aton/cn3  full  seal  e]  and  I 

A1  ( - )  {i  x  10^2  Btou/cH3  full  scale]  obtained  using 

the  progran  SLIC. 
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